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ABSTRACT 



Aims. How mass assembly occurs in galaxies and which process(es) contribute to such activity are some of the main questions highly 
debated in galaxy formation and evolution theories. This has motivated our survey MASSIV (Mass Assembly Survey with SINFONI 
in VVDS) of 0.9 < z < 1.9 star-forming galaxies selected from the purely flux-limited VVDS redshift survey. 

Methods. We evaluate the characteristic size and stellar mass of 45 MASSIV galaxies at 1 < z < 1.6 and we use the internal dynamics 
(~| , obtained with the SINFONI/VLT-ESO integral field spectrograph, in order to derive the stellar mass - size - velocity relations. We use 

^ ^ the Kennicutt-Schmidt formulation to estimate the gas content and to compute its contribution to the total baryonic mass in MASSIV 

I ■ galaxies. 

O [ Results. For the first time we obtain the relations between galaxy size, mass, and internal velocity, and the baryonic Tully-Fisher 

relation, from a statistically representative sample of star-forming galaxies at 1 < z < 1.6. We derive a dynamical mass which is in 
C/5 ■ good agreement with rotating galaxies containing a gas fraction of ~ 20%, that is perfectly consistent with the content derived using 

the Kennicutt-Schmidt formulation and corresponding to the expected evolution. Non-rotating galaxies are more compact in their 
extent of the stellar component, and less massive than rotators, but not statistically different in their gas extent. We obtain a marginal 
evolution in the size - stellar mass and size - velocity relations with discs being evenly smaller with cosmic time at fixed stellar 
mass or velocity, and less massive for a given velocity with respect to the local Universe. This result does not imply an abnormal 
evolution in the galactic spin as previously reported. The scatter of the Tully-Fisher relation is reduced introducing the 5 05 index, that 
' is interpreted with the increasing contribution to galactic kinematics of turbulent motions with cosmic time. We report a persisting 

scatter for rotators in our relations, that we suggest to be intrinsic, and possibly caused by complex physical mechanism(s) at work in 
our stellar mass/luminosity regime and redshift range. 

Conclusions. Our results consistently point towards a mild, net evolution of these relations, comparable to what is predicted by 
OA ' cosmological simulations of disc formation. In a conflictual picture where earlier studies reported discrepant results, our findings put 

, on firmer ground the lack of an influential transformation of the fundamental relations of star-forming galaxies for at least 8 Gyr and 

04 ■ a dark halo strongly coupled with galactic spectrophotometric properties. 



Key words. Galaxies: high-redshift - galaxies: evolution - galaxies: kinematics and dynamics - galaxies: fundamental parameters 



' 1. Introduction merging is the main physical process driving the assembly of 

. dark- matter halos and successfully describes the galaxy cluster- 

In the cuii-ent pai-adigm of structure formation, galaxies are properties, it is unclear how it affects the build-up of galax- 
formed m dark matter halos which grow along cosmic time ^^^^ ^nd how other physical processes might play a role, 
via hierarchical assembly of smaller units. While halo-halo 

With considerable efforts, large and deep spectroscopic sur- 



Sendojfprint requests to: D. Vergani, e-mail: vergani@iasfbo.inaf.it ^^^^ j^^^j higher-z galaxies have been performed in 

* This work IS based mainly on observations collected at the j^^^ ^^^^^^ understand galaxy formation and evolution. 

European Southern Observatory (ESO) Very Large Telescope, Paranal, . ^, ^j^u^ui.rji.^ 

Chile as part of the Programs 179.A-0823, 177 A-0837, 78.A-0177, "^^"^ f^""^', found between various 

75.A-0318, and 70.A-9007. This work also benefits from data products S^^^^y properties (e.g., luminosity, colour, stellar mass, surface 

produced at TERAPIX and the Canadian Astronomy Data Centre as brightness, star-formation rate) a general consensus has been 

part of the Canada-France-Hawaii Telescope Legacy Survey, a collabo- reached on a luminosity /mass-dependent evolution of the dis- 

rative project of NRC and CNRS. ti'ibu tion functions (e^. Zucca et aL 2006; Pozzetti et al. 20(37l 

** All the data published in this paper are publicly available following Bu ndv et al.1 120061: Montana et a l.1 |20M 120061) supporting the 

this link: http://cosmosdb.lambrate.inaf.it/VVDS-SINFONI, downsizing scenario (e.g. iCowie et al.lll996l) . 
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There is convincing evidence that 1 < z < 2.5 is a criti- 
cal cosmic time, when galaxy properties are evolving rapidly. 
During this period we witness an intense mass assembly activity 
along with a morphological differentiation of galaxies well un- 
derway. This epoch corresponds to the maxi mal star formation 
activity in the Uni verse around z ~ 2 (e.g. iTresse et aljF2007l; 
ICucciati et al.ll20fl ). and the fast build-up of the stell ar mass 
observed in early-type g alaxies at z ~ 1-1.5 ( Arn outs et alj 
120071: IVergarii et al.ll2008l) . It is also an epoch when the connec- 
tion between AGN and star formation could play an important 
role in shaping galaxies. In this framework the knowledge of the 
dynamical state of galaxies is an indispensable information to 
assess the mechanisms governing the growth of luminous struc- 
tures in the Universe. 

In the last years the availability of integral field spectro- 
graphs working under excellent seeing conditions, or assisted 
by adaptive optics, has opened a new era in the investiga- 
tion of galactic resolved dynamics at high redshifts, with in- 
tegral field spectrographs bringing considerably more informa- 
tion than traditional long-slit techniques. Mapping the full two- 
dimensional velocity fields and velocity dispersion is essential 
in many dynamical studies in contrast with the classical long- 
slit technique: it avoids misinterpretation of the measured ve- 
locities due to the misalignment between the photometric and 
kinematic axis, or to large dispersion velocities. This is espe- 
cially true at increasing look-back times when galaxies show a 
larger incidence of irregular shapes and chaotic motions. In ad- 
dition, new opportunities have recently been opened in the near- 
infrared regime in which new instruments allow one to sample 
the rest-frame optical emission that contains all major diagnos- 
tic lines, e.g., SIN FONI on VLT, Eisenhauer et al, (2003 ) and 
OSIRIS on Keck, iLarkin et alj (|2006|) . Spatially-resolved ob- 
servations of UV/optically selected galaxies at 1.5 < z < 3 
(Genzeletal. 2006, 2008; Forster Schreiber et al. 2006. 
Bouche et al. 2007; Law et al. 2007, 2009i IWright et all 



2009; 



2008 



Stark et al. 2008; Epinat et al. 2009) together with cosm ological 



simulations (e.g.. lGenel et al.1 120081; iDekel et alj|2009h indicate 
that many of the physical properties of these sources cannot be 
primarily attributed to major merger events experienced in the 
history of galaxies. Instead internal me chanisms are contributing 
to mass assembly dGenzel et alj|200§ , and references therein). 
In driving the observed dynamics, they can correspond to the 
infall of accreting matter, and/or large collisional clumps (e.g., 
iLaw et a l. 2007, 2009; Genzelet al. 2011) and/or processes re- 
lated with star formation (iLehnert et alJl2009l) . 

In summary the emerging scenario is minimizing the impor- 
tance of major merging in shaping up galaxies at early epochs, 
apparently at odds with the hierarchical assembly of mass in 
the cold dark matter paradigm - although merging systems are 
identified in about one third of the galax ies studied. A general 
conse nsus has not been reached yet (e.g.. lRobertson & BuUockl 
l2008h . 

So far a detailed view of dynamics of distant galaxies is re- 
stricted to a limited number of objects caused primarily by the 
large investment of telescope time to build up a fair sample. The 
most efficient approach to this issue is to adopt certain selec- 
tion techniques, such as various combinations of colours to fo- 
cus on a specific redshift range and/or targeting extended discs 
to ensure sufficient sp atial sam Eling (iLaw et alJ l2007l l2009t 
iForster Schreiber et alj 2006, 20(M). While pre-selecting galax- 
ies remains the most effective way to gather a large, high-z 
galaxy catalogue, when exploring statistical issues it is highly 
desirable to have a fully representative catalogue of galaxies with 
well known selection functions on multiple fields of a homoge- 



neous survey (e.g. JPuech et al.l2010h . Having a sufficient spatial 
resolution is also an issue as demonstrated by a survey such as 
IMAGES (Intermediate MAss Galaxy Evolution Sequence) that 
selects z ~ 0.6 galaxies observed with limited spatial resolu- 
tion (0'.'52 pix~^). The resulting picture emerging from various 
data-sets offers discrepant results on the evolution of the funda- 
mental relations between galaxy size, mass, and velocity (e.g., 
|Flores et al. 20 06; Puech et al. 2008, 2010; Bou che et al. 20o| 
Forster Schre iber et alj2009l;[Cresci et al. 2009; iLaw et alj2007l 
2009). As these quantities represent important indicators to un- 
derstand the structure and evolution of disc galaxies, the inter- 
pretation of their changes, or constancy in time, would imply 
a modification in the standard picture of galaxy disc formation. 
For instance, at odds with earlier findings and theoretical expec- 
tations, the absence of evolution in the size-velocity relation of 
disc galaxies reported by Bouche et al. (2007) have fundamental 
implications in the growth of the specific angular momentum of 
dark matter halos. 

To overcome the above-mentioned difficulties and shed light 
on relations linking various galactic properties to dynamics, we 
are carrying on the MASSIV programme (Mass Assembly Survey 
with SINFONI in VVDS) to obtain near-infrai-ed integral field 
spectroscopy of a hundred of star-forming galaxies selected with 
a pure flux limiting criterion on three VVDS fields. The main 
goal of MASSIV is to obtain a detailed description of galaxy dy- 
namics to probe their formation and evolution in a representative 
sample of galaxies in the crucial 1 < z < 2 epoch. In the first 
set of papers we present the survey in 
kinematic classification of galaxies in 



Contini et al"] (l201lh . flie 
Epinat etal.r ( l20 1 2l) . and 



their resolved metallicity in Oue vrel et al. ( 201 I h. In the present 



paper, we present the fundamental scaling relations that exist 
among the mass, size, and velocity in galaxies. 

In Sect. 2 we summarize our sample selection, observations 
and data reduction. In Sect. 3 we detail the procedures to de- 
termine the physical properties of the MASSIV galaxies (mor- 
phology, mass, star formation rate, kinematics inferred from the 
photometric and integral field spectroscopic data). In Sect. 4 we 
discuss the gas content and the size of our galaxies and in Sect. 5 
we summarize the impact of our resulting scaling relations in the 
context of mass assembly. 

Throughout this work we assume a standard cosmological 
model with Q.m = 0.3, Qa = 0.7 and Hq = 70kms 'Mpc ' 
(which gives a median scale of 8.32 kpc arcsec"' at z ~ 1.24). 
Magnitudes are given in the AB system. 



2. Observations 



2.1. Sample selection 



We report SINFONI dEisenhauer et al.f 2003VB onnet et al.ll2004l) 
observations of our new resolved kinematic MASSIV survey. The 
galaxies have been selected from the VVDS. The VVDS is a 
redshift survey which has mapped and presented results on the 
evolution of galaxies, large scale structures, and AGNs from the 
spectroscopic identification of about 50,000 source s down to a 
limiting magnitude of I^b - 22 in the Wide Fields (Garil li et al.l 
2008), to Iab = 24 in the Deep Field (Le Fevre et al. 2005|), and 
to Iab = 24.75 in the Ultra-Deep Field (see Tab.©- The VVDS is 
combined with a multi-wavelength dataset from radio to X-rays. 

The MASSIV sample analysed in this work was selected on 
the basis of several criteria (see details in Contini et al. 2011). 
They can be summarized as follows: 

1) Galaxies with a redshift assigned to their rest-frame UV 
VVDS spectra at a confidence level larger than 80% in the red- 
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shift range z = 0.9 - 1 .8 were pre-selected. The confidence level 
of the redshift determination has been estimated in the VVDS 
on the basis of repeated observations. Galaxies are grouped into 
four classes and we used only those that correspond to the max- 
imum success probability of having an accurate redshift. 

2) Their spectral features (Ha/ieSeSA, or [OlII].i5007A) had 
to fall in a region not contaminated by OH sky-lines. Strong sky- 
lines were outdistanced by a minimum of 9A from the above- 
mentioned spectral lines computed using our VIMOS spectro- 
scopic redshifts. Thus, we could observe an OH-uncontaminated 
emission line with a typical rotation of 200 km s"' . 

3) Star-forming galaxies below z < 1.46 have been se- 
lected in taking a certain threshold in the equivalent width (EW) 
of the [OlI]/l3727 doublet in VIMOS spectra with a signal-to- 
noise ratio (SNR) above 6. That is, z < 1.46 galaxies with 
EW[OlI] larger than 25A in SNR<10 spectra and 40A in 6 < 
SNR < 10 spectra have been selected. Star-forming galaxies 
at higher redshifts have been selected using the rest-frame UV 
spectra. These constraints guarantee the detection of HQ'/i6563 A 
(or [OlII]/i5007A) in SINFONI spectra within the scheduled ob- 
serving time. 

Out of the ~ 30, 000 spectra observed in the various VVDS 
fields, 10% of them satisfy these criteria (see Tab. [T] for the de- 
tails). From this subsample 84 galaxies have been observed in 
the framework of MASSIV. The galaxy VVDS220148046 sup- 
posed to be at redshift 1.3710 based on the VIMOS spectrum 
turns out to have a redshift of 2.2442 identified with the HyS 
and [OIII] lines in SINFONI observations. This galaxy is not in- 
cluded in the current analysis. Finally, our 83 galaxies are repre- 
sentative of the star-forming galaxy population at 0.9 < z < 1.8 
uniformly distributed in the stellar mass interval \ogMI Mq - 
9.47 - 11.77. In this sample ten galaxies have been observed 
with the adaptive optics (AO) system of SINFONI assisted with 
the Laser Guide Star (LGS) facility to achieve a better angular 
resolution. 

In this analysis we are using more than half of the en- 
tire MASSIV sample, i.e., 49/83 galaxies, or about ~ 60%. 
Thanks to the high precision of the VVDS redshifts, almost 
all observed galaxies have been detected (45/50) with either 
Ha^6563A (44/45) or [OIII]/15007A (1/45). In particular, with 
seeing-limited mode all but one (43 out of 44) galaxies tar- 
geted in Hq'/16563A have been detected. Seven out of 50 galax- 



Table 1. Sample selection. 



VVDS 


lAB.lim 


Number 


Criteria 


MASSIV This 


Survey fields 




of targets 


Full-filled 


sample work 


(1) 


(2) 


(3) 


(4) 


(5) (6) 


Wide 


17.5-22.50 


24,507 


224 (183) 


21 (4) 21 (4) 


Deep 


17.5-24.00 


12,668 2,600 (300) 


29 (1) 27 (1) 


Ultra-Deep 


23.0-24.75 


1,200 


86 (55) 


34 (6) 2 (2) 



(1) The fields of the VVDS used in MASSIV. (2) Ub limiting magni- 
tude selection. (3) Total number of spectroscopic targets. (4) Number 
of targets fulfilling criteria described in Sect. l2.1l (5) Total number of 
galaxies in the MASSIV sample at the completion of the survey. (6) 
Number of galaxies used in this analysis. 

In columns from (4) to (6) we indicate with parentheses the numbers of 
targets observed with adaptive optics in our observing campaign. 



ies have been observed with adaptive optics (out of which two 
have no detection, one targeted in Ha and the other one in 
[OIII]). The fifty galaxies that constitute the original MASSIV 
sample (see Tab. 1) have a median redshift of z = 1.2423 (where 
Z25% = 1 .0375 and Z75% - 1 .3547 are the two extreme quartiles 
of the redshift distribution). Excluding our serendipitous detec- 
tion at z = 2.2442 and the four undetected galaxies, the median 
redshift of our 45 galaxies analysed in this work remains simi- 
lar to the original sample (z - 1 .2246 with Z25% = 1 035 1 and 
Z75% = 1.3338). 

2.2. SINFONI Data Acquisition and Reduction 

The observations of the MASSIV sample have been collected 
with the near-infrared integral field spectrograph SINFONI 
of the VLT under the Large Programme 179.A-0823 (P.I. T. 
Contini) between April 2007 and January 2011 and including 
the observations of the pilot runs 75.A-0318 and 78.A-0177. Our 
subsample of 50 galaxies refers to observations conducted until 
2009 and fully reduced and analy zed before 201 0. 

In the accompanying paper bv lContini et alJ (|20I ih the strat- 
egy and global physical properties of MASSIV galaxies are de- 
tailed, along with a more exhaustive explanation of selection 
fu nctions. Full details on the data reduction procedures are given 
in lEpinat et alj ( |20I2|) . Here we provide a short summary. 

We used the J- and H-bands to sample the spectral interval 
1 .08 - 1 .41 fim and 1 .43 - 1 .86 fim with a spectral resolution of R 
~ 2000 and ~ 3000, respectively. The majority of our observa- 
tions were in seeing-limited mode with a 8"x 8" field-of-view 
and 0'.'125 x 0'.'25 pixel scale. For a subsample of seven galaxies 
we obtained AO observations with a 3'.'2 x 3'.'2 field-of-view and 
O'/IOO X 0'.'050 pixel scale. To maximize the telescope time, we 
offsetted the majority of our targets at different positions in the 
field-of-view of the instrument. This strategy avoids sky frame 
acquisition (see for details in Contini et al. 201 1). In addition, we 
applied a sub-dithering to place the target at a different position 
on the chip. Conditions were photometric during the observa- 
tions with a median seeing of 0'.'70 measured on the PSF stars 
acquired during the night. To accurately point our galaxies, we 
acquired them through a blind offset from a bright nearby star. 
We also observed standard stars over the same nights to flux cal- 
ibrate. Individual exposures range between 300s and 900s with 
a total on-source integration time that ranges between Ih to 2h. 

The core of the data reduction has been perf ormed using the 
ESO-SINFONI pipehne, version 2.0.0 (Modi ghaniet al.ll2007l) 
complemented with additional routines to homogenize the data 
processing among the different reducers. We subtract the sky 
background, skylines, and dark cuiTent from the raw science 
data using the contiguous sky frames obtained with the offset-to- 
sky sequence. The data were flat-fielded using an internal lamp, 
wavelength calibrated from the spectra of reference arc-lamps. 
These processed data frames were then reconstructed into data 
cubes. Individual cubes of a given observing block were aligned 
in the spatial direction by relying on the telescope offsets from 
a nearby bright star and then combined. The rejection of cosmic 
rays has been applied on the combined cubes. The observations 
of a telluric standard star of a typically B spectral type follow the 
science frames and were reduced in a similar way. We extracted 
the mono-dimensional spectra of the stars by summing the flux 
of an aperture of diameter of 5 resolution elements to obtain the 
flux calibration after removing the telluric absorption lines. We 
coiTected for the atmospheric transmission dividing the science 
cubes by the integrated spectrum of the telluric standard. We 
created sky cubes to quantify the effective resolution for each 
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science data cube. The sky cubes were reduced in a similar way 
as science frames, but no correction for sky subtraction has been 
applied. We fitted with a Gaussian profile the unblended night 
sky lines on the extracted, mono-dimensional spectra. The effec- 
tive FWHM spectral resolution has a typical value of 130 km s"' . 



3. Derived quantities 

3.1. Morphology and kinematics 

We use Ha^6563A (and [OIII]/15007A Une in one case) to de- 
rive the dynamical properties. We assume that the ionized gas 
rotates in a thin disc with a velocity reaching a plateau in the 
outer regions (but for 8/23 rotating galaxies marked with ^'^^ in 
column (9) of Tab. 2 this radius is not reached within the area 
covered by our observations). Using ax^ minimization we pro- 
duce seeing-corrected velocity and dispersion maps of galaxies 
with geometrical inputs weighted by the SNR of each pixel. 

The input geometrical parameters to the fitting model are es- 
timated from the I-band best-seeing CFHTLSQ Megacam images 
for all galaxies, but 14^' Wide galaxies that were covered with 
the I-band CFHT-12K camera only (McCracken et al. 2003). At 
the typical z ~ 1.2 redshift of MASSIV galaxies, I-band images 
probe the 3200-4200A rest-fra me wavelength ra nge, or U-band. 
We use the GALFIT software (iPeng et al.l 2002h that convolves 
a PSF with a model galaxy image based on the initial fitting 
parameter estimates with a Sersic (1968) profile. GALFIT con- 
verges into a final set of parameters such as the centre, the posi- 
tion angle, and the axial ratio. Residual maps from the fit were 
used to optimize the results. The I-band images were also used 
to calibrate for SINFONI astrometry, in considering the relative 
position of the PSF star. GALFIT has also been applied on the 
ionised gas maps (properly deconvolved with a PSF star) in order 
to derive the semi-major axis disc scale-length. The CFHT im- 
ages and the kinematic maps are presented in Epinat et al. (2012) 
where more discussion is presented on the modeling procedure. 

The major source of uncertainties in deriving the maximum 
rotation velocity is the inclination of the studied galaxy. For very 
small objects the inclination is not robustly constrained, thus for 
a fraction of galaxies the median inclination value of a randomly 
distribution of galaxies on the sky (60°) is adopted instead. For 
this category of objects we have a typical uncertainty of 24° as- 
sociated to the inclination. Another uncertainty appears when the 
maximum velocity is not reached within the radius covered by 
our observations. Other uncertainties entering in the error bud- 
get are related to the simplistic disc model adopted. To evalu- 
ate the net effects of this assumption would require the assump- 
tion of different, detailed modeling schemes, which is beyond 
the scope of this paper The final uncertainty associated to our 
velocities includes both the errors originating from the deprojec- 
tion of the radial velocity given a certain inclination (estimated 
with Monte Carlo simulations) and those from the modeling pro- 
cedure (qua ntified with simulations on the GHASP local sample 
foUowing lEpinat et al.ll2012h . 



' Based on observations obtained with MegaPrime/MegaCam, a joint 
project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii 
Telescope (CFHT) which is operated by the National Research Council 
(NRC) of Canada, the Institut National des Science de TUnivers of the 
Centre National de la Recherche Scientifique (CNRS) of France, and 
the University of Hawaii. This work is based in part on data products 
produced at TERAPIX and the Canadian Astronomy Data Centre as 
part of the Canada-France-Hawaii Telescope Legacy Survey, a collab- 
orative project of NRC and CNRS. 



3.2. Mass and star formation rate 

In this work we obtain the stellar mass from spectral energy dis- 
tribution (SED) fit to photometric and spectroscopic data with 
BC03 stellar population synthesis models (Bruzual & Chariot 
200 31) using the GO SSIP Spectral Energy Distribution tool 
(iFranzetti et al.ll2008l) . We assume a Salpeter initial mass func- 
tion (Salpeter: 1955ir with a lower and upper mass cutoff of re- 
spectively 0.1 and IOOAIq, a set of delayed exponential star 
formation histories with galaxy ages in the range from 0.1 to 
15 Gyr. As input for the SED fitting we use the multi-band 
photometric observations available in the VVDS fields, includ- 
ing BVRI data from the CFHT-12K camera, ugriz data from 
the CFHT Legacy Survey, J and Ks-band data from SOFI at 
the NTT, from the UKIDSS survey, SW IRE data when avail- 
able, and the VVDS spectra. Following IWalcher et aP (l2008h 
we adopt the probability distribution function to obtain the stel- 
lar mass (listed in column (6) of Tab. 2), the absolute magni- 
tude, and the other results of fitting procedure as detailed in 
Contini et al.. (.201 1) . 

We estimate the dynamical mass and the gas mass in our 
galaxies to trace the evolution of both baryonic and dark mat- 
ter components. As described in the next section our galaxies 
show a variety of kinematics that we classify in two broad cate- 
gories: rotators and non-rotating system s (see details in Sect. 3.3 
of this paper and in lEpinat et al.ll20T2l) . The dynamical mass is 
estimated using the following equation: 



Mdyn 



y2 R 0-2 
~G~ Gh^ 



(1) 



where the first term represents the mass enclosed within the 
radius R defined as the total extent of the ionized gas, V is the 
velocity at the same radius, and G is the universal gravitational 
constant. The second term represents the asymmetric drift cor- 
rection following lMeurer et al.l (119961) where cr is the dispersion 
velocity at R, h is the gas disc scale length in an isotropic ge- 
ometrical configuration with a surface density described by a 
Gaussian function. The dynamical mass of each galaxy is given 
in column (6) of Tab. 2, its error is estimated accounting for the 
uncertainties related to the inclination and the modeUng proce- 
dure. 

The gas mass is derived with the empirical correlation be- 
tween star formation rate (SFR) and gas surface density assum- 
ing that the gas and the stellar content accounts for the majority 
of the mass in the inner centre of galaxies. In the local Universe, 
iKennicut? ('1998b') proved that the surface density of gas to star- 
formation rate follows a power law above some critical gas den- 
sity (Schmidt 1959, I,sfr = ^^"as) over more than six orders of 
magnitude in Ssfr- This empirical correlation has not yet been 
routinely proved at high redshift because of the lack of exten- 
sive measurements of cold gas mass, but a number of studies 
on small sam ples support the validity of this relation at high 
redshi fts (e.g., 'Erb et al."2006'; 'Bouche et al. 2()07t iDaddi et alJ 
120101: iTacconi et al.n2010) . Some detailed studies on individual 
galaxies also point to a consistency with the local Schmidt law 
(e.g., the lensedz = 2.7 LBG MS1512-cB58. lBaker et all2004 . 
Thus we compute the star formation rate surface densities using 
the total SFR and the size of our galaxy. 

The SFRs are computed following iKennicuti (Il998al) . The 
fluxes are measured on > 2cr regions of the flux-calibrated data 
cubes and corrected for dust reddening using the extinction coef- 
ficient derived from the SED fitting (see Contini et al. 201 1). The 
half-light radius, R\/2,gas, has been me asured as the sem i-major 
axis half-light radius with GALFIT dPeng et all l2002h on the 
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Table 2. Properties of the MASSIV sample 



ID 

(1) 



z 

(2) 



Line 
(3) 



M, 



(4) 



gas 

(5) 



(6) 



Rlasl 

(7) 



^I/2,sni 

(8) 



R 



1 /2,.«ar 

(9) 



10.09!«;| 

10.10!; 

10.18!| 

10 13+''" 

10.08!«]^ 

10.72^11 

10 15+"^" 

10.17+"!^ 

10.02+1 ° 

10 02+" 
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0.72+ 0.09 
1.31+0.25 
2.68+ 0.21 
3.96+ 0.12 
2.66+ 0.11 
4.04+ 0.16 
1.93+ 0.19 
1.80+ 0.40 
1.39+ 0.98 
5.45= 
8.96= 
2.62= 
2.07= 
2.89= 
3.32= 
2.76+ 0.25 
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Tab. 2. (1) VVDS identification number. (2) Redshift derived from SINFONI observations of the emission line (Ha.l65 63A or [OIII]/1500 7A) 
tabulated in column (3). (4) Stellar mass in solar mass units. (5) Gas mass derived using the Kennicutt-Schmidt formulation dKennicutll I998bl) . (6) 
Dynamical mass computed following Eq. 1 and based on the best-fitting parameters of the kinematic model. (7) The radius computed as the total 
size of the ionized emission map with a confidence level larger than 2cr. (8) The gas half-light radius estimated on the SINFONI maps of the lines 
tabulated in column (3). (9) The observed I-band half-light radius. (10) Maximum velocity deduced from model fitting. (11) Ratio of the maximum 
rotation velocity over the velocity dispersion. (12) Kinematic classification (R: Rotating objects, NON-R: non-rotating systems). " This redshift is 
computed on the original VIMOS spectrum. '' Observations with Adaptive Optics. '■' Galaxy with low signal-to-noise ratio (3 < SNR < 4.5) plotted 
with empty symbols in the figures. Galactic properties of SINFONI observations with a signal-to-noise ratio below 3 are not reported. '' Galaxy 
for which the plateau velocity has not been reached. ° Galaxy for which the criterion on v/cr > 1 is not fulfilled, thus not included in the purely 
rotating disc sample. The public access to data is at ,http://cosmosdb.iasf-milano.inaf.it/VVDS-SINFONI^ 
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spatial extent of the ionized gas (Ha/ieSeSA or [OIII]/15007A) 
emission, after convolving with the PSF of a bright star nearest 
in time to the scientific observations. Maps of the residuals were 
used to optimize the results. Assuming the Schmidt law holds 
in our redshift interval we estimate the gas mass Mgas using the 
relation between the gas density Sgas an d the star-formation rate 
surface density Esfr (iKennicutnl 1 998bh : 



-SFR 



M0 yr^' kpc"2 



= 2.5x10" 



-gas 



Mopc- 



1.4 



(2) 



or 



Non-Rotating systems: At least one of the conditions on 
the residuals, position angles, or rotation-over-chaotic mo- 
tions is violated. In this class fall galaxies with some kind of 
anomalies in the velocity field, merger-like systems, and/or 
objects with off-peaked distributions of the dispersion ve- 
locity, but also possibly slow rotators or face-on systems 
with the observed velocity gradient lower than typically 
+25kms"'. 



where the relation between the observed Ho- luminosity per 
unit area, or the Ha surface density Sh», and the gas surface 
density is: 



M0 pc- 



= 1.6 X 10 



-27 



x0.71 



,erg s ' kpc 



(3) 



The mass of the gas associated with the measured star for- 
mation is given by: 



Mgas = 2^gas X Area 



(4) 



where the Area in Eq. 4 is the total spatial extent of the ion- 
ized gas measured on > 2cr regions of the flux-calibrated data 
cubes. 

The gaseous masses are given in Tab. 2, the associated error 
is analytically propagated from the errors in the observations us- 
ing conventional techniques. Further details on the error analy - 
sis are described in the accompanied paper dEpinat et al.ll2012h . 
Typical value of the error is of the order of 0.25 dex. We note 
that diff'erent calibratio ns exist for the Kenn icutt-Schmidt law 
at higher redshift, e.g., iBouche et al.l (l2007h suggest a steeper 
power-law index (1.7 instead of 1 .4 in Eq. 3) for z > 2.5 starburst 
galaxies. This steeper formulation implies a < 0.3 dex difference 
in the gas computation. Given the uncertainties involved in de- 
riving this quantity, a different value for the slope does not have 
any significant impact on our data interpretations (see Gnerucci 
etal. 2011). 

3.3. Kinematic classification 

To classify our sample we adopt a kinematic classification that 
re lies on a combinatio n of observed galaxy properties. We refer 
to lEpinat et al.l (1201 2h for a full description of the methodology 
for which we provide here a brief outline. In a first step eight 
of us classified independently all the sample using personal cri- 
teria. Subsequently we reconciled the criteria into a new agreed 
scheme: we assign to each galaxy a confidence flag generated by 
the concordance among the classifiers. 

When galaxies satisfy all the following conditions they are 
kinematically flagged as: 

Rotating discs: when the velocity field is well described by a 
symmetrically rotating disc with a kinematic position angle 
that does not differ significantly from the morphologically 
derived major axis. We quantify these conditions assuming 
that the mismatch between the position angles is lower than 
+20° and that the residuals of the modelled velocity field 
over the measured velocities are less than a certain fraction 
(< 20%). In this work we also require that galaxies classified 
as "rotators" are rotationally supported imposing that the ro- 
tation velocity must be larger than the dispersion velocity 
(i.e., v/cr > 1 as listed in column (12) of Tab. 2.) 



4. The gas properties in massiv star-forming 
galaxies 

We compare the total extent of the ionized gas, Riasi, with the 
half-light radius of the stellar component, Ri/2,siar, for the class 
of rotating and non-rotating disc galaxies in Fig.[I] (in the top 
panel with blue colour-coded circles and in the bottom panel 
with red colour-coded squares, respectively). Riasi is the radius 
computed as the total size of the ionized emission map that has a 
confidence level larger than 2cr. The half-light radius of stars is 
measured on the observed I-band images, or rest-frame U-band 
wavelengths. The error bars in Fig.[T]are the Icr uncertainties de- 
rived from the GALFIT procedure on the images corrected for 
their respective PSF FWHM. The histograms show the distribu- 
tions of the gas and stellar radii of the classes. The solid line is 
the linear fio to the data and the dashed line shows the standard 
deviation. 

The dotted line is the correlation observed in a thin expo- 
nential disc between the radius of the total ionized gas emis- 
sion and the half-light radius of the stellar component, Riast ~ 
1.9 X R\/2,siar- Assuming that the optical radius equals 3.2 times 
the disc scalelength following Persic & Salucci ( 1991), the half- 
light radius is 1 .678 times the disc scale length. The total ionized 
gas, Riast, is c lose to the optical r adius as in local star-forming 
galaxies (e.g. JGarrido et al.ll2005h . 

For the class of rotators we obtain median values of Riast = 
5.404gy kpc and R\i2,star - 3.5 1 2 62 kpc, and for non-rotators 
R,as, = 5.50^:^° kpc and Ryi2,star = 2.68^:^^ kpc (with the ex- 
tremes representing the quartiles of the distribution for the two 
classes). Non-rotating galaxies are more compact in the extent of 
the stellar component than rotators, but similar in the gas extent. 
For rotating galaxies we obtain that = 1.63 + 0. 14x7?i/2,.5/flr 
with an intrinsic scatter of cr,„,r = 2.23 (and total scatter of 
cr,ot - 2.70). The same correlation for non-rotating systems is 
Riast - 1-70+0.21 x/?i/2,jrflr(o"m/;- - 2.20 and cr,o; - 2.66). In the 
local Universe the extent of the ionized gas com ponent in star- 
forming galaxies is close to the optical radius (e.g.. lGarrido et al] 
2005). Within Icr dispersion, and considering the associated 
uncertainties, all MASSIV rotating galaxi es are within the ex- 
pected statistical trend found bv Persic & Saluccil(ll991h . A sim- 
ilar agreement has been found by Puech et al. (2010) for their 
IMAGES sa mple at z - 0.6. Similarly to Puech et al. (2010) 
and Bamford et al.l (l2007l) . we notice a fraction of non-rotating 
galaxies with a larger gaseous extent than their stellar counter- 
part exceeding the predicted correlation by Icr. While in the 
sample of Puech et al. (2010) the cases where this correlation 



- We fit a linear relationsliip to tiie data using thie MPFITEXY rou- 
tines (Williams etal. 2010). The MPFITEXY routine depends on the 
MPFIT package (Markwardt 2009.) . This routine properly adjusts the 
intrinsic scatter to ensure the minimisation of chi-square with an itera- 
tive prescription. 
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is not respected are the ones with the UV light extending fur- 
ther out the field-of-view of the FLAMES/IFU (3" x 2"), in 
the MASSIV sample the two outliers show a signature of a close 
companion, one of which is detected in Ho'. 

We will further investigate these trends and their relation 
with the environment at the completion of the survey when we 
will develop a classification scheme enabling us to distinguish 
the various mechanisms of disturbance in inclined and face- 
on galaxies as well as in compact and dispersion-dominated 
spheroids. For simplicity, at the present time we flag all these 
galaxies in a unique class named non-rotators (plotted with red- 
coded, square symbols in all figures). 

We now compare our measurements with the characteristic 
sizes of galaxies as found by other similar surveys to analyse the 
correlation among properties like stellar mass, size, and veloc- 
ity. In IMAGES (Puech et al. 2010) the total Ught radius mea- 
sured on the ionized [Oil] gas of rotating galaxies has a median 
of 11.29 kpc and ranges from 8.25 kpc (first quartile, or 25th 
percentile of distribution) to 11.99 kpc (last quartile, or 75th 
percentile). The median value of the half-light rest-frame UV 
radius is 6.27 kpc with 4.26 kpc and 6.53 kpc as quartiles of 
the distribution. We note that the IMAGES disc size is the max- 
imum between the total size measured on the ionized [Oil] gas 
and the half-light rest-frame UV light multiplied by 1 .9. In SINS 
(Cresci et al. 2009), the half- width half-maximum size (HWHM) 
measured on the ionized Ha maps is interpreted as exponential 
disc scale-lengths of z ~ 2 rotating SINS discs, and ranges from 
4. 1 1 kpc to 6.35 kpc with median of 5.87 kpc. The circular half- 
lig ht sizes of the total SIN S z > 2 sample computed on Ha maps 
bv'Forster Sc hreiber et all (12009) is i n the range 2.4 — 4.60 kpc 
and has a median value of 3.10 kpc. Dutton et alJ (1201 li) claim 
for a fundamental disagreement of the SINS HWHM sizes de- 
rived by Cresci et al. ( 2009) and the SINS half-light radii by 
iForster Schreiber et alJ (l2009h . as already noted by these latter 
authors who reject the possibility to attribute the discrepancy to 
non-e xponential discs. At hi gher redshift the AMAZE/LSD sam- 
ple of lGnerucci et alJ (|201 ih is constituted by eleven z ~ 3 rotat- 
ing galaxies. The characteristic radius of these exponential discs 
measured on Ha maps ranges between 0.71 kpc and 1.73 kpc 
(first and last quartile of the distribution) with a median value 
of 1 .24 kpc. Overall, comparing directly to the values of Puech 
et al. (2010) as we adopt the same radius definition, we find 
that IMAGES radii are a factor of 1.5 larger than MASSIV radii. 
Instead, compared to those of SINS by Forster Schreiber et al] 
(l2009l) . we find similar sizes (our half-light sizes measured on 
He maps are in the interval Ri/2.gas = 2.85 - 3.96 kpc with 
a median of 3.47 kpc). Both SINS and MASSIV radii are ap- 
proximately three times larger than sizes sampled at z ~ 3 in 
AMAZE/LSD by Gnerucci et al. (201 1). 

We investigate further the disc properties to verify the hy- 
pothesis whether our MASSIV discs have exponential profiles. 
Discs of ionized gas have frequently an irregular, clumpy dis- 
tribution as it is also reported for SINS galaxies. Despite the 
simplistic assumption both the Sersic index distribution and the 
correlation between exponential and half-light radii suggest that 
MASSIV discs can be approximated to exponential discs. The 
same qualitative conclusion has been reported for SINS discs 
(iForster Schreiber et alJl2009l) . 

In summary, our resulting characteristic size agrees within 
the statistical correlations previously assessed in the literature. 
The radius taken to derive the velocity is critical to avoid the 
insertion of systematic biases (see Noordermeer & VerheiienI 
l2007h . The relation obtained between the size of the stars and 
that of the gaseous component ensures that this latter physical 




3 6 9 12 

Rl/2,,tar [kpc] 



Fig. 1. Comparison between Riast and Ri/2,star radii computed as 
the size of our ionized gas maps with a confidence level larger 
than 2cr and the half-light radius of the stellar continuum mea- 
sured on the observed I-band best-seeing CFHTLS images, re- 
spectively. Within the errors both rotators (top panel, blue cir- 
cle) and non-rotators (bottom panel, red square) agree with the 
expected statistical trend Riast = 1.9 x R\i2,star found by Persic 
& Salucci (1991, PS91) plotted with a black dotted Hne. The 
blue/red solid lines are the best-fit to the rotators/non-rotating 
galaxies (dashed lines show the total scatter of the correlations). 
The empty, dashed symbols show galaxies detected with a low 
signal-to-noise ratio (3 < SNR < 4.5). 



quantity is reliable for the determination of the maximum veloc- 
ity. 

Figure|2] shows the stellar mass of MASSIV galaxies versus 
the dynamical (top panel) and gas (bottom panel) mass. The 
methods used to derive these masses are described in Sect. l3.2l 
The total distributions of these quantities are shown as full gray- 
coded histograms while those with blue, positive-oblique angle 
(or red negative-oblique angle) refer to rotating (non-rotating) 
galaxies. Dynamical masses range between 2.57 x 10'° Mo and 
1.29 X 10" Me with a median value of 5.01 x 10'" Mq for the 
class of rotating galaxies (dashed-point horizontal line in the bot- 
tom panel of Fig.©, and 0.95 - 8.13 x 10'" Mq with a median 
value of 2.45 x 10'" Mo for non-rotating galaxies (dashed hori- 
zontal line). We estimate a gaseous mass in the interval between 
3.38 X 10'' Me (first quartile) and 1.64 x 10'° Mo (last quartile) 
with a median value of 7.85 x 10^ Mo, and stellar masses in 
the following range 1.48 - 5.87 x 10'" Me with a median value 
of 3.98 X 10'" Mo in rotating galaxies. Median values for non- 
rotators are M,,,,^ = 1.23 x 10'° Me and M,,,, = 6.71 x 10'' Me 
with ranges between AM,,^, = 0.90 - 1.90 x 10'" Mo and 
AM^„j = 3.51 -9.37 X 10'' Mo- 

Within their uncertainties only few galaxies lie in the forbid- 
den region of M^tar > Mdyn- The only strongly deviating rotating 
galaxy (empty, dashed symbol) has been detected with a low 
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Fig. 2. The stellar mass content of MASSIV galaxies compared 
to the dynamical mass derived from the dynamical modeling 
(top panel) and to the gas content derived using the Kennicutt- 
Schmidt formulation (bottom panel). Symbols are as in Fig. 1. 
The horizontal and vertical histograms represent the distribution 
of the stellar and dynamical masses (top) and gas mass (bottom 
panel). The horizontal and vertical lines are the median values 
of rotators (point-dashed blue line) and non-rotators (dashed red 
line), respectively. Dashed diagonal lines in the bottom panel 
show the condition: Mgas - frac x Mstar where frac is equal to 
25% for rotating galaxies and 42% for non-rotators (or 20% and 
30% of the baryonic ma ss, respectively). The solid line is the lo- 
cal relation proposed by[Schi^novic3(200^ using the GALEX 
Arecibo SDSS survey of cold HI gas. In the inserted panels 
the cumulative distributions for stellar-, dynamical-, and gas- 
mass are plotted for rotators (solid, blue line) and non-rotators 
(dashed, red line) along with the probability for the two classes 
to be statistically different. 



signal-to-noise (~ 3.5). We have a probability of f = 7.5 x 10"^ 
and P - 1.5 X 10 "* that the stellar mass is uncorrected to dy- 
namical and gas mass using the Spearman correlation test. The 
correlation is statistically significant in particular for the stellar 
and gas masses. In the inserted panels of Fig.|2]we show the cu- 
mulative distributions of these quantities for rotators (solid, blue 
line) and non-rotators (dashed, red line) along with the proba- 
bility for the two classes to be statistically different using the 
Kolmogorov-Smirnov test. 

We assume no contribution from dark matter in the inner 
centre of galaxies. It is the best hypothesis to assume given the 
inability to measure the dark matter distribution in high-z galax- 
ies. This hypothesis has been adopted in earlier studies, e.g., by 
Gnerucci et al. (2011). With this assumption the mass in gas 
taken as difference between dynamical and stellar mass of both 
classes of rotator and non-rotator galaxies, is on average a frac- 
tion of ~ 45% of the dynamical mass (diagonal dotted line), 
and approximately 20% of the stellar mass or 17% of the bary- 
onic mass. A consistency picture is obtained for rotating galaxies 
when deriving the gas mass from the Kennicutt-Schmidt law. We 
fit a content of gas that is 25% of the stellar mass (or 20% of the 
baryonic mass) in rotating galaxies with typical errors associated 
to the fit of the order of 5 - 6%. 

The gas content in non-rotating galaxies is definitively 
higher (~ 42% of the stellar mass, or 30% of the baryonic mass) 
when using the Kennicutt-Schmidt law. Taking into account that 
this class of galaxies may contain face-on discs, this fraction 
represents a lower limit. The same value derived as difference 
between dynamical and stellar mass is somehow lower. It can 
be justified by two factors: 1) the inclusion of merging systems 
in the class of non-rotating galaxies that may have not yet ex- 
hausted their gas reservoir; 2) the different properties of dark 
matter halos in rotators and spheroidal galaxies. In nearby galax- 
ies, within the optical radius the dark matter content is suspected 
to be higher in rotators than in spheroids (Persic et al. 199<^. 
Thus for the same dynamical mass, the fraction of gas and stars 
should be lower in rotators than in our class of non-rotating sys- 
tems. 

If we fix the s lope to the local relation proposed by 
ISchiminovichI (1200 8') using the GALEX Arecibo SDSS survey 
of cold HI gas, we obtain a weak evolution of +0.17 dex to 
z = 1.2 for the entire sample (or -i-O.ll dex for rotators and 
-1-0.21 dex for non-rotating galaxies). We emphasize that the se- 
lection criteria adopted to build up the MASSIV sample lead to a 
subsample of non-rotating galaxies with a lower content of stel- 
lar mass at fixed gas and dynamical mass. 

The fraction of gas in z ~ 0.6 IMAGES galaxies ranges be- 
tween 30% as derived from the evolution of the gas-metallicity 
relation (R odrigues et al. 2008) and 45% of the stellar mass ob- 
tained with the inverse Kennicutt-Schmidt law by Puech et al. 
(2010). The gas fraction quoted by Cresci et al. (2009) for SINS 
galaxies at z ~ 2 is in the interval 23 - 30% of the dynami- 
cal mass, but including a dark matter contribution of 40%. At 
higher redshift (z > 3, the AMAZE/LS D surveys) these values 
are too uncertain dGnerucci et al ]l201ll) . 

Based on direc t estimates using mo lecular observations (e.g., 
iDaddi et al ] 120 lot iTacconi et al.1 12010) . the gas mass in disc 
galaxies ranging in the interval between 34-44% of the bary- 
onic mass, is slowly decreasing in the last 8-10 Gyr Our results 
are in full agreement with a slow decrease of the gas mass con- 
tent in disc galaxies since z ~ 1 .2 with a gas fraction that halves 
progressively down to the local Universe. 
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5. The size - mass - velocity relation 

In this section we examine the fundamental relations existing in 
the MASSIV sample at 1.0 < z < 1.6. We compare their sizes, 
stellar masses, and rotation velocities to galaxy samples avail- 
able in the literature at different redshifts to investigate the evo- 
lution of these properties. In particular, we derive the relation 
between the luminosi ty (or stellar mass) and rotational veloc- 
ity as introduced by ITullv & Fisheii (11977b . or Tully-Fisher re- 
lation (TFR). We further explore the shift in the zero point of 
this relation and the other fundamental relations that unveil the 
inter-connection between dark and luminous matter for galaxies 
at different cosmic epochs. 



5.1. The size - mass relation 

The evolution of the size - stellar mass relation in rotating galax- 
ies at 1 .0 < z < 1 .6 is shown in the top panel of Fig.|3] The distri- 
bution of the MASSIV rotating galaxies in this plot (blue-coded 
circles) is consistent with a mild dependence between the disc 
size and the stellar mass. Taking into account the relatively lim- 
ited statistics associated to a disc size representing only a lower 
limit (as measured on galaxies with a signal-to-noise ratio above 
3), we interpret this plot in the frame of its evolution to the lo- 
cal relation fixing the slope to the local relation as derived by 
[Outtonet al. (2011). This relation is based on the local I-band 
disc scalelength - stellar m ass relation from Dut ton et al. (2007) 
with s tellar masses from iBell et al.l (l2003h and iBell & de Jong! 
(l200l . The evolution of the relation is only slightly stronger us- 
ing the V-band instead of the I-band size-stellar mass relation 
derived from the SDSS data. This radius is scaled to the half- 
light radius which is 1.678 times the scale length for the expo- 
nential profile. This local relation is plotted in Fig. 3 as dashed, 
black line. The solid, blue line is the fit to the MASSIV rotating 
galaxies. To make a consistent comparison to our measurements 
we plot in this figure only the evolution reported in those studies 
using the sizes measured on the ionized gas (He, or [Oil]) maps. 
The best fit obtained using SINS galactic sizes defined as the 
HWHM Ha radii (Bouche et al. 2007; Cresci et al. 2009) is plot- 
ted with a dotted-dashed, orange line and with long-dashed, or- 
ange line when the size is defi ned as the half-light Ha radius by 
iForster Schreiber et alJ (l2009h . Note that difference in the evo- 
lution of this relation is due to the method used to measure the 
disc sizes by [B ouche et al. (2007); Cresci et al. (2009) and by 
iForster Schreiber et alJ (120091) (cf. lDutton et al.ll201 li) . We plot 
with a dotted, red line the evolution measured on IMAGES ro- 
tating galaxies using as disc size the maximum between the total 
size measured on the ionized [Oil] gas and the half-light rest- 
frame UV light multiplied by 1.9. All quantities are consistently 
rescaled to our Salpeter IMF. 

There is a general consensus on the mild evolution of 
the disc size in star-fo rming galaxies up to z ~ 1 (for ex- 
ample as measured by W illiams et al.l (|2010) taking the cir- 
cularized rest-fra me I-band ha l f-light radii of disc-dominated 
galaxies, and by iTruiillo et al.l (l2006l) taking the circularized 
rest-frame V-band half-light radius). Compared to the present- 
day disc sizes, the z = 1 counterparts are smaller by ~ 
0.10 dex at fixed stellar mass and rotation velocity (e.g., 
iDutton et al.l 1201 1\ Th is picture agrees well with the theoret- 



ical predictio ns ( Blum enthal et al. 19841; Dalcanton et al. 19971; 

' 1 nn o. n-; „™ — : p. A,.:i„ r> ^"onnnL lj„ t p_ t J 



Mo etal."l998; Firmani & Avila-Reese"200(y; I de Jong & LacevI 
2000; Pizagno et al. 2005; Dutton et al. 2007). 

At z ~ 1.2 we report an evolution size-stellar mass relation 
of -0.14 dex using the half-light Ha radius that is quite simi- 



lar to earlier findings at z ~ 1. Above z ~ 1 the evolution of 
disc sizes in literature is less clear, even with opposite trends. 
Trujillo et al. (2006) found a factor of two smaller half-light 
sizes at z - 2.5 at fixed stellar mass relative to z = galaxies. 
A similar evolution is reported by Williams et al. (2010) in star- 
forming galaxies at z = 2. The evolution of cir c ular h alf-light 
Ha sizes as measured bv Forste r Schreiber et al.l (l2009h and re- 
ported by Dutton et al. (201 1) on SINS data at similar redshifts 
is smaller than previously quoted results, i.e., -0.07 ± 0.05 dex 
to be compared with -0.28 dex by Trujillo et al. (2006). This 
discrepancy is attributed by Dutton et al. (2011) to the differ- 
ence in the half-light radius measured on the Ha maps and rest- 
frame optical band images. Taking into account this difference 
their results can be reconciled. A different trend is obtained us- 
ing measurements of Cresci et al. (2009) on the same SINS data 
by Dutton et al. (201 1). Using the half-width half-maximum size 
interpreted as exponential disc scale-lengths, SINS disc galaxies 
presented in Cresci et al. (2009) are ~ 1.6 times larger at z = 2 
than their local counterparts at fixed stellar mass. Even if a fac- 
tor of -0.07 dex caused by the skewed distribution towards low 
inclinations reduces the discrepancy, the values remain signifi- 
cantly inconsistent by a factor of three with the measurements 
of Trujillo et al. (2006) and WilHams et al. (2010) (cf. Dutton et 
al. 201 1). This discrepancy between data-sets increases at higher 
redshifts if we explore the AMAZE/LSD sample by Gnerucci et 
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Fig. 3. The stellar mass-size relation (top panel) and the 
velocity-size relation (bottom panel) for MASSIV rotating 
galaxies at z ~ 1.2 (solid, blue line with the Icr error of the 
correlation illustrated by the light blue area). The local relation 
by Dutton et al. (2007, D07) is shown as short-d ashed, black 
line. W e overplot the best fit relations obt ained by Puech et aP 
( 120081) (P08. z ~ 0.6, do tted red line), bv lBouche e t al. (2007h 
and ICresci et all (l2009l) at z ~ 2 (B07, dotted-dashed or- 
ange line), and by Dut t on et al. (2011) using measurements by 
IForster Schreiber et alj(l2009h (FS09, z ~ 2, long-dashed orange 
line). The symbols are as in Fig. 1. 
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al. (201 1). The sizes of the discs are more than four times higher 
in the past than in their local counterparts. We can cautiously 
speculate that this strong evolution may be caused by a selection 
effect in the AMAZE/LSD galaxies and the difficulties in mea- 
suring realistically the Ha sizes at such high redshifts (z ~ 3). 

Comparing the evolution computed directly from the half- 
light radii measured in Ha, we obtain at z ~ 1.2 the similar 
evolution as derived from SINS data a t z ~ 2 by Dutton et 
al. (2 011), based on measurements by iForster Schreiber et alj 
( l2009i) . However if we account for the shift suggested by Dutton 
et al. (201 1) that converts our Ha sizes at z ~ 1.2 to their optical 
band half-light radii, we should have a stronger evolution with 
redshift, i.e., -0.21 dex instead of -0.09 + 0.04. Conversely, we 
do not find the same conversion as suggested by Dutton et al. 
(2011) in the ratio of galaxy sizes (e.g., between our observed 
1-band and Ha). Furthermore, our results on the evolution of 
the disc size is perfect ly in ag reement with the evolutio n up t o 
z ~ 1 found bv .WilUamset all f2010) and lTruiillo eTal] (l2006h . 
especially taking into account the associated uncertainties on the 
measurements and the adopted assumptions. Therefore based on 
the current data-set we suggest a mild evolution in the size - 
stellar mass relation of rotating galaxies at z ~ 1.2. An homoge- 
neous comparison of MASSIV galaxies at z ~ 2 and z ~ 1 that 
will be possible at the completion of this survey, will provide a 
definitive answer to this controversial issue. 

5.2. The size - velocity relation 

The bottom panel of Fig. 3 shows the relation between the disc 
scale length and the maximum rotation velocity for MASSIV ro- 
tating galaxies and for samples from the literature at different 
redshifts. 

The distribution of the MASSIV rotators in this plot (blue- 
coded circles) is consistent with no dependence between the 
disc size and the stellar mass. For the reasons explained in 
Sect. 5.1, we discuss hereafter the size-velocity relation using 
the slope determined for local galaxies. We observe a shift of 
-0.12 + 0.05 dex of the relation at z ~ 1.2 with respec t to the 
local 1 -band size-veloci ty distributio n computed by Dutto n et al.l 
(l2007h . As reported bv lPuech et alj 12007) different calibrations 
proposed to describe this correlation show very similar results, 
despite the different quantities used to trace the galactic proper- 
ties. Puech et al. (2007) found some IMAGES rotating galaxies 
at z ~ 0.6 with a relatively lower disc scale length at fixed ve- 
locity compared to local rotators (cf. blue circles in their Fig. 3), 
although the majority of them lies on the distribution of local ro- 
tating discs. We fit their rotating galaxies with a relation showing 
a marginal evolution of -0.08 ± 0.06 from z ~ 0.6 to z ~ 0. 

At z ~ 2 the results differ with authors. A consistent evolu- 
tion with our measurements is observed at higher redshift adopt- 
ing th e half-light sizes as computed by ' Forster Schreiber et alJ 
(l2009l) . The evolution at z ~ 2 shows evenly smaller sizes at 
a given maximum rotation velocity relative to lower redshift 
galaxies with an offset of -0.22+0.06 dex from z = 2 to z = (cf . 
Dutton etal. 201 1), and -0.10 dex to z ~ 1.2. At z ~ 2adifrerent 
trend is observed by Bouche et al. (2007) combining data from 
the SINS survey (Cresci et al. 2009) and from Courteau ( 1997) 
at z = 0. Bouche et al. (2007) found no evolution of the zero- 
point of the size-velocity relation from z = 2 to z = 0. Dutton et 
al. (2011) attribute this discrepant result to an inconsistency in 
the SINS half-width half-maximum size reported in Cresci et al. 
(2009) and interpreted as exponential disc scale-lengths. 

The size-velocity relation shows at all cosmic times a large 
scatter (of the order of ~ 0.2 dex). The scatter reported in previ- 




ses [kms-' ] 

Fig. 4. (Top) The stellar mass TF (smTF) relation at z ~ 1.2 
based on the MASSIV rotating galaxies. Symbols are as in Fig. 1 . 
Errors on the velocities are computed using GHASP simula- 
tions to account for various uncertainties. The smTF for the 
MASSIV rotators is calibrated using either the local slope defined 
by Pizagno et al. (2007, P07) or by Bell & de Jong (2001, BdJOl) 
(dashed black line and solid black line, respectively). The solid- 
blue lines are the best-fit to MASSIV rotating galaxies and the 
cyan shaded area shows the Icr error on the zero-point param- 
eter (the intrinsic and total scatter is plotted as cyan dotted line 
and cyan dashed line, respectively). Other fits are by Puech et aP 
(20Q§) (P08, z ~ 0.6, dotted red fine), Cresci et al. (2009 ) (C09 
z ~ 2 .2, short-dashed-dotted orange line), and iGnerucci et al.l 
('2011') (Gll,z ~ 3, long-dashed-dotted magenta line). (Bottom) 
The relation between the stellar mass and the contribution of 
both ordered and chaotic motions (5 05 = ^/{Q.5 x v^„, + cr^), 
Kassin et al. 2001, KOI) is shown. The cyan shaded area shows 
the Icr error on the zero-point parameter, the cyan dotted line and 
cyan dashed line are the intrinsic and total scatter, respectively. 
We overplot the best-fit relation by Kassin et al. (2007) at z ~ 1 
(KOI, dot-dashed black line) and bv iGnerucci et al.l (1201 II) at 
z ~ 3 (Gl 1, dot-dashed magenta line). Symbols are as in Fig. 1. 
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ous analyses at high-z are likely due to the difficulties in discern- 
ing kinematically complex systems, that translate in uncertain- 
ties in the estimation of the characteristic radius of the system. 
While we are not immune to this subtle possibility, the statistical 
trend observed between the size of our ionized emission maps 
and our radius of the stellar continuum on the broad-band im- 
ages of rotating galaxies, ensures reliable quantitie s. The local 
relations studied by e.g., Courteau ( 1997) and Mat hewson et alj 
119921) and better constrained in that respect, also show a large 
scatter of this relation suggesting a physical meaning of this dis- 
persion. 

5.3. The mass - velocity relation 

In the top panel of Fig. 4 we show the correlation between the 
stellar mass and the rotation velocity at z ~ 1 .2 of the MASSIV 
galaxies. Rotating galaxies are plotted with blue-coded circles. 
They are the sole class of galaxies included in our best fitting 
procedure. We plot non-rotating galaxies with red-coded, square 
symbols. The errors on the velocities are computed using sim- 
ulations on a local sample to account for various uncertainties 
(see for details Epinat et al. 2008). The errors on the stellar mass 
are described in iContini et ah (i201 1) . Consistent to what has 
been done in the previous section, given the relatively limited 
statistics at this stage of the survey we set the slope to the lo- 
cal calibration. To make a fair comparison with the published 
data at different redshifts, and to allow for future comparisons, 
we fit our data with the two most widely used calibrations. In 
particular, we calibrate the stellar m ass TF relation (smTFR) 
with the local slope bv iPizagno et al . (2007) to compa re wit h 
the IMAGES sample at z ~ 0.6 bv jPuech et al.. (.2008^ l2010h . 
The local calibration bv lPizagno et al.l (l2007h based on a repre- 
sentative galaxy subsample extracted from the SDSS and revised 
by [Hammer et al. (2007 ) is plott ed in Fig. |4] with dashed, black 
Une. The best fit bv lPuechet al.l(l2008l) is plotted with a dotted 
re d line. We calibrate t he relation with the local slope defined 
bv lBeU & de Jond (1200 1) to compare with the SINS sample at 
Z ~ 2.2by|Cresci et al. (2009) and with the LSD/ AMAZE galax- 
ies at z ~ 3 b y Gnerucci et al. (2011). The local calibration by 
iBell & de Joii g (2001 ) is plotted i n Fig.Hwith a solid, black line. 
The best fit bV iCresci et al.*('2009) is s hown with a s hort-dashed- 
dotted orange line and that by Gneruc ci et al.l(l201 ih with a long- 
dashed-dotted magenta line. The Bell & de Jong (2001) calibra- 
tion is based on the local sample by Verheiien (2001) who se- 
lected predominantly l ow-mass, gas-rich galaxies observed wit h 
the 21cm HI emission dPuech et al.ll2010tlHammer et al.ll2007l) . 
The resulting relation is steeper than Pizagno et al. (2007) refer- 
ence. The two relations start to be equivalent around a velocity 
of ~ 200 kms"' and a stellar mass of ~ 10" Al©, thus samples 
with these galactic properties are not substantially affected by 
the choice of the calibration. The vast majority of our MASSIV 
galaxies are located below these values. In consequence, adopt- 
ing one, or the other calibration aff'ects the strength of the evolu- 
tion in the velocity-mass relation. We compute for each calibra- 
tion the zero-point of the correlation, its error, its intrinsic and 
total scatter All the above-mentioned studies have been oppor- 
tunely rescaled to account for our Salpeter IMF. 

For rotators we obtain an evolution from z = in stellar mass 
of the smTFR zero point of -0.36 ± 0. 1 1 dex, using as reference 
the SDSS subsample bv iPizagno et al.l (l2007l) (with an intrinsic 
scatter of cr,-,,,^ = 0.32 dex and a total scatter of cr,,,, - 0.48 dex). 
Instead, using the Bell & de Jong (2001) calibration we obtain a 
zero-point at z ~ 1.2 that is consistent with no evolution since 
z ~ (or -0.05 + 0.16 dex, cr,„,, = 0.52 and o-,o, = 0.72). 



The solid, blue lines are the best-fit to MASSIV rotators and the 
cyan shaded area shows the 1 cr error on the zero-point parameter 
(the intrinsic and total scatter is plotted as cyan dotted and cyan 
dashed line, respectively). 

The shift in the zero point of the smTFR with respect to z = 
computed bv lPuech et a l. (2008, 2010) at z ~ 0.6 is surprisingly 
high (-0.34 dex in stellar mass) if compared to our results at 
z ~ 1.2. The evolution in the smTFR computed using IMAGES 
dataset is very similar (within the errors) to the evolution we re- 
port at higher redshift (-0.36 ± 0.1 1 dex at z ~ 1.2). There are 
important issues to be taken into account w hen interpreting th ese 
values. The original evolution claimed by Flores et aP (l2006h on 
the same data-set was consistent with no, or marginal evolu- 
tion. Other authors, using traditional slit spectroscopy, reached 
the same qualitative conclusion of no/mild evolution up to red- 
shift z - 1 (e.g., Weineretal. 2006; Conselice et al. 2005|; 
[Fernandez Lorenzo et al.l 120101) . The IMAGES galaxies have 
large seein g and coarse spatia l resolution due to the FLAMES 
instrument ( Puech et alJl2010h . Furthermore, the lack of evolu- 
tion in the baryonic (star-i-gas) Tully-Fisher relation reported 
by Puech etal. (2010) using the same IMAGES sample and a 
gas fraction at z ~ 0.6 similar to the values of nearby galaxies, 
would imp ly a minimal, if not negligible evolution of the stellar- 
mass TFR dDutton et al.ll201 11) . All these points would therefore 
suggest an over-estimation of the evolution of the stellar Tully- 
Fisher relation by Puech et al. (2008) and an evolution of the 
smTF relation of IMAGES data of the order of -0.1 dex only, 
instead of the estimated -0.34 dex (iDutton et al.ll201lh . 

If we compare our best fit with the Bell & de Jong (2001) 
calibration at z ~ 1.2 of the four SINS galaxies by Cresci et 
al. (2009) at similar redshifts, it is consistent with no evolu- 
tion of the zero point. These four galaxies are very massive 
(~ 10" Al©) with velocities around ~ 250 km s"' where the two 
calibrations give similar estimates. At larger redshift, Cresci et 
al. (2009) found a shift in the stellar masses at fixed velocity 
of -0.41 +0.11 dex compared to the z = Bell & de Jong 
(2001) calibration. Similar evolution is observed with the cal- 
ibration by Pizagno et al. (2007), or -0.44 ±0.11 dex (Cresci 
priv. comm.). The absence of discrepancy when taking two dif- 
ferent calibrations is expected giving the properties of the anal- 
ysed SINS galaxies. They are very massive and fast rotators and 
are located at the intersection between the two local relations in 
the stellar mass - velocity diagram. Gnerucci et al. (201 1) using 
LSD/AMAZE galaxies show the possible build-up of the stel- 
lar mass TF relation at z ~ 3. Despite the large scatter of the 
correlation (~ 1.5 dex) they report a statistical shift of the zero 
point to their local relation of - 1 .29 dex (a very fast evolution, 
-0.88 dex, between z ~ 2.2 and 3 compared to more recent cos- 
mic epochs). 

Galaxies at high redshifts have a larger veloci ty dispersion 
and a higher fraction of turbulent motions (e.g., iGenzel et al.l 
2008). To account for the large velocity dispersion that con- 
tributes to galactic kinematics at increasing redshift, we plot in 
the bottom panel of Fig. |4] the relation between the stellar mass 
and the 5 05 index. The 5 05 index is defined as follows: Sqs — 
V(0.5 X v^^, + cTq), it accou nts for the contribution of both or- 
dered and turbulent motions dKassin et al.l2007tlCovington et al.l 
I2OIOI) . 

The cyan shaded area shows the Icr error on the zero-point 
parameter, the cyan dotted and cyan dashed lines are the intrinsic 
and total scatters, respectively. Diagonal lines represent the best- 
fit relation by Kassin et al. (2007) at z ~ 1 (dot-dashed black 
line) and by .Gnerucci et al.i (i201 li) at z ~ 3 (dot-dashed ma- 
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genta line) calibrated using the slope bv lKassin et alJ (l2007h . We 
observe a tighter correlation with a reduced scatter than the clas- 
sical smTF relation. Physically this smaller scatter of the distri- 
bution reinforces the recent results on the importance of chaotic 
motions to the galactic kinematics, and the larger contamination 
encountered in traditional long-slit spectroscopy. 

Using the 5 05 index the rotators and non-rotators are lo- 
cated now in the same locus (bottom panel of Fig. 4). The 
z ~ 1.2 smTFs()5 is log(M,^)= 4.46±0.07 + 2.92 x log(So5) 
(with cr,„„.=0.08 and a-,oi-QA5) when galaxies of both classes 
are included. If only rotating galaxies are considered (this is 
the correlation plotted in the figure), the relation can be writ- 
ten as log(M*)= 4.26+0.10 + 2.92 x log(So5) (with cr,„,r=0.13 
and cr,o,=0.43). We therefore report a shift in stellar mass 
of -0.25 dex of the rotating MASSIV sample (o r -0.45 dex 
includ ing both dynamical classes) compared to iKassin et alj 
(200"/) g alaxie s at z ~ 1 and of -1-0.35 dex (or +0.55 dex) to 
lOnerucci et al. galaxies at z ~ 3. 

Figure 5 shows the relation between the baryonic mass and 
the velocities (top panel) and the ordered+chaotic motion com- 
ponents, quantified using the 5 0.5 index (bottom panel). Based 
on our results we infer no evolution; that is to be compared 
with the conclusions reached by McGaugh] (l2005i) in the local 
Universe implying a very marginal evolution of the content of 
the gaseous component with respect to nearby galaxies. Previous 
results on the baryonic TF relation using a data-set of integral 
field survey have been derived by Puech at al. (2010) who find no 
significant evolution compared to the local Universe. Based on 
this issue these authors claim for gas already in place at z ~ 0.6 
with no need to advocate external gas accretion being the gas al- 
ready bound to the gravity well of galaxies. Other interpretations 
are possible: if the growth with cosmic time of galactic veloci- 
ties and masses happen along the scaling relation, as predicted 
by LCDM models, any apparent evolution can be observation- 
ally recognized. 

6. Discussion and Conclusions 

We have explored the relationships between galaxy size, mass, 
and internal velocity of the statistically representative MASSIV 
sample of 45 galaxies with near-infrared resolved kinematics at 
1 < z < I.6., a n unbiased sample of high-z star-forming galaxies 
dContini et al.„ 2011). 

Before the advent of the 3D integral field spectroscopy 
the possible sources of discordance in using these quan- 
tities include observational diffic ulties and selection effects 
dPortinari & Sommer-Larsen '^2007). The higher fraction of mor- 
phological and kinematic disturbances in high-z galaxies com- 
pared to nearby objects may bias the real evolutionary effects. 
However, systematics and subtle selection effects also need to 
be taken into account when comparing with recent near-infrared 
resolved kinematics surveys. For example, the observations of 
the OSIRIS sample (iLaw et al.ll2007l 120091 IWright et alj|2008l) 
were assisted by adaptive optics to reach very high angular res- 
olutions, but probing only the highest surface brightness emis- 
sion. With this observing mode they sample smaller characteris- 
tic Wa radii (~ 1.3 kpc) than those typically probed (~ 3 kpc) in 
othe r similar survey s at c omparable re dshifts. The IMAGES sur- 
vey dHammer et al.il2007t|Pue"ch et alj|2008ll2010l) samples both 
star-forming and more quiescent galaxies, but the coarse spatial 
resolution of FLAMES instrument requires meticulous correc- 
tion methods to extract the maximum velocity. The SINS sam- 
ple is the largest collection of galaxies at high-z, but extracted 
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Fig. 5. The baryonic TF (bTF) relation versus the maximum ve- 
locity (top panel) and versus the v + cr (or 5 0.5) index (bottom 
panel). We overplot the best-fit relation of rotating galaxies (blue 
symbols) using the local slopes by McGaugh (2005) (plotted as 
solid black line in the top panel, McG05) and by Kassi n et"al] 
(l2007h at z ~ 1 (plotted as solid black line i n the bottom panel , 
K07). The other fit in the top panel is by iPuech et alJ d2010h 
(PIO, z ~ 0.6, sohd, orange line) and in the bottom panel by 
i Gnerucci et al.l d201 ih (z ~ 3, long-dashed-dotted magenta line, 
Gl 1). The cyan shaded area shows the Icr error on the zero-point 
parameter (the intrinsic and total scatter is plotted as cyan dot- 
ted line and cyan dashed line, respectively). Symbols are as in 
Fig.l. 
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Fig. 6. Star Formation Rate versus stellar mass for rotators (left 
panel) and for non rotators (right panel) in the MASSIV (1 < z < 
1.6, blue circle), IMAGES {z ~ 0.6, red square), SINS (z ~ 2.2, 
orange triangle) and LSD/ AMAZE samples (z ~ 3, magenta as- 
terisk). The grey dotted lines show the empirical main sequence 
of star-forming galaxies following Bouche et al. (2010) at dif- 
ferent redshifts. All quantities have been opportunely rescaled 
to account for our Salpeter IMF. 



from surveys with different selection criteria. The SINS sample 
used to analyse the relationship between stellar mass and veloc- 
ity, represents galaxies detected with the highest signal-to-noise 
ratio: they are the most luminous/massive, fastest, and highly 
star-forming galaxies. 

How some of these observational facts reflect on the physi- 
cal interpretation is illustrated in the left panel of Fig. 6 where 
the star formation rate of rotating discs obtained by SED fit- 
ting in recent resolved kinematics surveys is plotted as a func- 
tion of the stellar mass. The class of non-rotating/merging sys- 
tems is shown in the right panel of the same figur e. The main 
sequen ce of star-forming galaxies as defined by B ouche et al.l 
(I2OIOI) is drawn at diff'erent redshifts. The IMAGES rotators lie 
below the sequence defined by the MASSIV rotators consistently 
rescaled in star formation rate and stellar mass given their later 
cosmic epoch (but covering a large interval of galactic proper- 
ties), similarly high-z LSD/ AMAZE rotating discs stay above the 
sequence. SINS rotators are highly star-forming objects and very 
massive, thus coincide with the massive, star-forming tail of the 
MASSIV galaxy distribution. The SINS properties can explain 
the mostly absent scatter in the Tully-Fisher relation by Cresci 
et al. (2009) compared to our scatter in the same relation (cf. our 
Fig. 4 with their Fig. 5). Contrary to what was expected, the dif- 
ference in SFR between MASSIV and SINS galaxies cannot be 
attributed to the higher redshift and/or to t he limits imposed o n 
the minimum line flux to ensure detection (IContini et al.ll20lTI) . 
In this figure there is a larger dispersion of SINS rotators around 
the main sequence compared to other samples. This fact may 
be the consequence of at least two combined effects. The first 



one may be due to the heterogeneous SINS selection. The sec- 
ond effect may be physical, caused by an unexpected larger in- 
cidence of internal/external mechanisms (inflows/mergers, out- 
flows/feedback) acting on both the spectrophotometric and dy- 
namical properties of rotating discs. While we cannot discard 
the second hypothesis, it is unlikely that the same mechanisms 
acting on SINS galaxies at z ~ 2 are not at play a few Gyr later in 
MASSIV galaxies (assuming both samples are representative of 
the entire rotating disc population). Nevertheless, we can spec- 
ulate that we witness the same phenomena acting with differ- 
ent efficiency on the various physical regimes (of galactic stellar 
mass, star-formation rate, etc.). Furthermore, comparing the SFR 
for rotators (left panel in Fig. 6) and non-rotators (right panel) 
for SINS and MASSIV surveys, it appears that non-rotating sys- 
tems display similar ranges of SFR at z ~ 1 and z ~ 2, while 
it is not the case for rotating systems that exhibit higher SFR at 
z ~ 2 than at z ~ 1 . Our results highlight the importance of the 
disc formation in the SFR mechanics between z ~ 2 and z ~ 1 
(assuming these samples are representative of the entire rotating 
disc population). In summary, data-sets available so far (even 
those taking advantages from resolved spectroscopy), are spec- 
troscopic surveys that pass through a natural pre-selection from 
one, or several photometric catalogues. As a consequence the 
comparison between samples remains difficult. At present, firm 
conclusions based on comparisons with other surveys, that are 
differently plagued by systematics in the galaxy properties de- 
terminations (e.g., different selection effects, quantities obtained 
at different wavelengths and with different assumptions), should 
be formulated with caution. 

In this work we aim at exploring statistical issues. We anal- 
ysed a large representative catalogue of galaxies with well- 
known selection functions, that is based on flux-limited multi- 
ple fields of a homogeneous redshift survey observed at higher 
spectral resolution, but with sufficient spatial resolution probing 
typical scales of ~ 5 kpc. Our results are the following. 

The MASSIV galaxies show a correlation between dynami- 
cal and stellar mass with an offset interpreted as a fraction of gas 
mass of the order of 20-25% of the stellar mass at z ~ 1 .2. A sim- 
ilar amount of gas mass is obtained using the Kennicutt-Schmidt 
formulation. This fraction of gas mass evolves mildly for the 
last 8 Gyr It is reinforced by the lack of a statistically significant 
evolution of our baryonic Tully-Fisher relation. Instead, the dis- 
agreement in the gas estimates using the two approaches for non- 
rotating galaxies, suggests different properties of dark matter ha- 
los in rotators and systems not supported by rotation at high-z, 
similarly to what is prescribed for nearby objects. The stellar 
mass in MASSIV galaxies is correlated with the incidence to be 
a rotator, or a non-rotator, with non-rotating galaxies presenting 
a lower content of stellar mass at fixed gas content and dynami- 
cal mass, but not at fixed star-formation rate. Furthermore, non- 
rotators are more compact in the extent of the stellar component 
than rotators, but not in the gas extent. 

We obtained correlations between galaxy size, mass, and in- 
ternal velocity, as observed in the local Universe and indicated 
by previous results, but the size - mass and size - velocity rela- 
tions are weak at 1 < z < 1.6. While individual galactic proper- 
ties are in fact evolving with time as advocated by most recent 
results of large galaxy surveys, the marginal change observed in 
the relationships linking the three fundamental quantities anal- 
ysed in this work - size, mass, and internal velocity - sug- 
gests an evolution along these relations. This picture agrees well 
with cosmological N-body/hydrody nam ical simulations by e.g., 
iPortinari & Sommer-LarsenI (l2007h and iFrrmani & Avila-Reesd 
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(l2009l) . We do not confirm the hypothesis of a strong, positive 
evolution in the size - stellar mass and size - velocity relations 
with discs being evenly smaller with look-back time at fixed stel- 
lar mass or velocity. Our results do not imply an abnorma l evolu- 
tion in the galactic spin as previously reported by B ouche et alj 
(l2007h . However, we discover a large spread in the distributions, 
the scatter is quite reduced for non-rotating galaxies when in- 
cluding the contribution of turbulent motions (with the 505 in- 
dex) in the relations involving the velocities. It is in disagree- 
ment with the hypothesis that large scatters in the low-velocity 
regime are due t o distu rbed or compact objects, as suggested by 
ICovington et al.l ([2010). In fact, restricting our analysis to reg- 
ularly rotating discs that sample a large range of stellar masses 
and sizes, we report a persisting scatter. Its origin is likely intrin- 
sic, and possibly caused by complex physical mechanism(s) at 
work in our stellar mass/luminosity regime and redshift range. 
At the present stage of the MASSIV survey and given the typi- 
cal scale (~ 5 kpc) used to trace the resolved kinematics, we are 
not able to investigate further the mechanism(s) driving these 
galactic properties. It is still a matter of debate whether we trace 
more turbul ent discs where star formation is t he energetic driver 
(e.g., lLehnertetalJl2009t iGreen et alJ lIoToh . or we are influ- 
enced by the infall of accreting m atter, a nd/or large co Uisional 
clumps (e.g.. iLaw et"al]|2007l 1200 9; Genz elet al.ll2011l) . and/or 
merger events (^g.- lHammer et al.ll201 10 . 

Further progress will be possible at the completion of our 
survey, doubling the statistics with galaxies in our highest (z > 
1.2) redshift range. However, some issues will require larger- 
scale surveys of kinematically resolved galaxies. Future surveys 
should be able to assess firm conclusions about the cosmic evolu- 
tion of the slope of the size/mass/velocity fundamental relations 
of discs at z > 1 . There is also the impelling request to establish 
the imprint of the environmental effect at work, preferentially at 
certain stellar mass regimes. Of course, the stellar component 
and ionized gas content of galaxies offer a partial view of the 
galaxy formation and evolution, which must be supplemented 
with the properties of the cold/molecular phase of the interstel- 
lar component. We need to comprehend under which conditions, 
and how efficient, the gas is converted into stars through cosmic 
epochs. The relation between the galactic properties probed in 
our survey and the characteristics of the gas reservoir will offer 
a crucial step to link the data to models of galaxy formation. 
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